Sublethal effects of parasitism on ruminants can have cascading consequences for ecosystems

Significance We found that pervasive parasitic infections reduce herbivory rates and can trigger trophic cascades. Lethal parasites clearly have cascading impacts on ecosystems, but whether common sublethal infections have similar effects is largely unknown. Using a mathematical model, we probed how parasites that reduce host survival, fecundity, or feeding rates can indirectly alter producer biomass in a helminth–ruminant system. We found that both lethal and sublethal infections triggered trophic cascades by altering the biomass of ruminant herbivore hosts and their resources. However, a global meta-analysis revealed that helminths tend to have pervasive sublethal effects on free-living ruminants, including by reducing host feeding rates. Our findings suggest there are widespread, but overlooked, ecological consequences of sublethal infections in natural ecosystems.

Parasitic infections are common, but how they shape ecosystem-level processes is understudied. Using a mathematical model and meta-analysis, we explored the potential for helminth parasites to trigger trophic cascades through lethal and sublethal effects imposed on herbivorous ruminant hosts after infection. First, using the model, we linked negative effects of parasitic infection on host survival, fecundity, and feeding rate to host and producer biomass. Our model, parameterized with data from a welldocumented producer-caribou-helminth system, reveals that even moderate impacts of parasites on host survival, fecundity, or feeding rate can have cascading effects on ruminant host and producer biomass. Second, using meta-analysis, we investigated the links between helminth infections and traits of free-living ruminant hosts in nature. We found that helminth infections tend to exert negative but sublethal effects on ruminant hosts. Specifically, infection significantly reduces host feeding rates, body mass, and body condition but has weak and highly variable effects on survival and fecundity. Together, these findings suggest that while helminth parasites can trigger trophic cascades through multiple mechanisms, overlooked sublethal effects on nonreproductive traits likely dominate their impacts on ecosystems. In particular, by reducing ruminant herbivory, pervasive helminth infections may contribute to a greener world.
helminth j herbivore j trophic cascade j feeding rate j producer Predators are well known to trigger trophic cascades by affecting the densities and traits of their prey (1). Emerging evidence suggests that parasites, including microparasites (bacteria, fungi, protozoa, and viruses) and macroparasites (helminths and arthropods) also trigger trophic cascades in similar ways (2,3). Both predators and parasites can have cascading effects on ecosystems by killing their victims [i.e., consumptive effects on density (2,4)] or by causing sublethal harm to victims prior to exploitation [i.e., nonconsumptive effects on traits (5,6)]. However, unlike predators, parasites also modify victim traits by consuming but not killing their hosts (2). This raises the question of whether consumptive (i.e., postinfection) impacts of parasites occur primarily through lethal or sublethal effects on hosts (2,7). Since nearly all organisms have longlasting, intimate associations with a multitude of parasites, understanding the extent to which parasites trigger trophic cascades may be equally, if not more, important than for predators (8). A better understanding of the effects of parasites is especially critical because changing environmental conditions are currently affecting the abundance and diversity of parasitic species across the globe (9)(10)(11)(12).
It is no surprise that lethal infections can have cascading effects on lower trophic levels and ecosystem functioning. One of the best examples of this was revealed by the eradication of rinderpest virus. In the late 19th century, rinderpest killed 80 to 90% of domestic and wild ruminants in sub-Saharan Africa. Following a successful vaccination campaign in cattle, release of wildebeest populations from infection dramatically altered carbon storage in the Serengeti ecosystem (13). However, unlike the rinderpest example, many infections are not lethal. Indeed, most animals survive countless infections caused by a range of parasites over their lifetimes. Such sublethal infections may also trigger trophic cascades through the modification of host behaviors or other phenotypic traits (7). For example, parasitic infection is associated with shifts in host body size and life history (14), movement behavior (15,16), dietary preferences (17), nutrient processing (18,19), vulnerability to predation (20)(21)(22)(23)(24), offspring sex and size (25), and both increased (26) and decreased (27)(28)(29) feeding rates. Changes in any of these host traits are likely to affect lower trophic levels, particularly if a host is abundant or plays a keystone role in the community.
Given the wide-ranging consumptive effects of parasites on hosts, parasite-initiated trophic cascades could emerge through a variety of mechanisms. Yet, the potential for parasites to indirectly affect ecosystem-level processes such as biomass production or

Significance
We found that pervasive parasitic infections reduce herbivory rates and can trigger trophic cascades. Lethal parasites clearly have cascading impacts on ecosystems, but whether common sublethal infections have similar effects is largely unknown. Using a mathematical model, we probed how parasites that reduce host survival, fecundity, or feeding rates can indirectly alter producer biomass in a helminth-ruminant system. We found that both lethal and sublethal infections triggered trophic cascades by altering the biomass of ruminant herbivore hosts and their resources. However, a global meta-analysis revealed that helminths tend to have pervasive sublethal effects on free-living ruminants, including by reducing host feeding rates.
Our findings suggest there are widespread, but overlooked, ecological consequences of sublethal infections in natural ecosystems. nutrient cycling is largely missing from eco-epidemiological theory (2), with notable exceptions for primary producer and bacterial hosts (30)(31)(32). This is because epidemiological models typically treat hosts and parasites as separate from the surrounding ecosystem, thereby ignoring other important components of the food web, such as predators and resources (33)(34)(35). Existing resource-host-parasite models indicate that lethal parasites can stabilize resource-host oscillations, regulate host density, and relieve herbivory pressure on plants, the latter of which causes trophic cascades (36). Infections that alter the relative biomass at different trophic levels (15,27) may also affect overall ecosystem stability (30,36). Although sublethal infections are incredibly common (37), the potential for parasiteinduced changes in host traits to have cascading consequences for ecosystems has rarely been explored theoretically (see ref. 29 for theory centered around microparasite infection in an aquatic invertebrate host). Likewise, there are very few empirical studies linking sublethal effects of parasites to trophic cascades, and those from terrestrial systems are almost exclusively limited to invertebrate hosts (refs. 2, 7, and 38 and references therein). It is therefore unclear if such findings would generalize to vertebrate hosts, their associated parasites, and terrestrial ecosystem-level impacts.
Helminths are important parasites of vertebrates, and the individual-and population-level effects of helminths on their hosts have been studied across a range of taxa (e.g., refs. [39][40][41][42][43][44]. In this study, we focus on free-living ruminants (order Artiodactyla), because they are a globally distributed group of mammals that occur in almost every terrestrial biome. These hosts occupy central positions in terrestrial food webs as key primary consumers and prey for secondary consumers, and many have large impacts on ecosystems through herbivory and nutrient cycling via waste (45). Because of the economic importance of both domesticated and free-living ruminants, their parasites are well-described (46). Likewise, the consumptive impacts of parasites, including helminths, on ruminants have been intensively studied (47,48). Importantly, infection by helminths can be lethal (e.g., refs. 39, 49, and 50) or sublethal, and the latter can have both reproductive and nonreproductive consequences for ruminant hosts (e.g., refs. 25, 41, and 51-53). Together, these traits make ruminants and their helminth parasites an excellent system to investigate how lethal vs. sublethal (including reproductive and nonreproductive) effects of infection scale up to affect host population dynamics and ecosystem function.
In our study, we combined mechanistic models with a comprehensive meta-analysis to investigate the potential for parasitic infections to trigger trophic cascades via their effects on a suite of key host traits. First, using mechanistic models, we tested the extent to which helminth-induced variation in the survival, fecundity, and feeding rate of herbivorous hosts affects the population dynamics and biomass of parasites, hosts, and primary producers. To explore how parasite-induced changes in these host traits might trigger trophic cascades, we parameterized our model using published data from caribou and reindeer (both Rangifer tarandus, henceforth referred to as caribou) and their common gastrointestinal helminth parasites, Ostertagia spp. (48). Caribou and their helminths are among the beststudied wild ruminant-parasite systems, in part because of the ecological, economic, and cultural importance of caribou in tundra ecosystems (54)(55)(56)(57) and the ongoing disease threats to the system posed by climate change (48,58). Second, we conducted an independent and comprehensive meta-analysis of the effects of helminth infection on the same traits (survival, fecundity, and feeding rate) of free-living ruminant hosts in order to evaluate empirical support for the mechanisms by which helminths are most likely to trigger trophic cascades in natural ecosystems globally. We also quantified the effects of helminth infection on host body mass and body condition, because changes in these traits could impact host survival, reproduction, and feeding rates (e.g., refs. 59-61). Our findings reveal that helminth parasitism has the potential to trigger trophic cascades via both lethal and sublethal effects on hosts but that sublethal effects are dominant in nature. As a consequence, sublethal parasitic infections may be an overlooked driver of ecosystem-level processes.

Results
Model Results. To explore effects of infection in a macroparasiteherbivorous host system, we used a model that builds upon a previous framework by Smith and Grenfell (36) to track densities of producers (A), herbivorous hosts (H), adult parasites within hosts (P), and free-living stages of the parasite in the environment (E) through time using ordinary differential equations ( Fig. 1). We used estimates from the literature for the genus Ostertagia (Table 1) to parameterize the model for caribou hosts that are commonly infected by Ostertagia gruehneri and that experience negative effects of these infections (51,62). Caribou diets consist of a mix of plants and lichens. In our model, we represent these autotrophic resources with a single state variable because they are both being consumed; for simplicity, we refer to them throughout as producers. First, we used a fixed set of parameters to evaluate the stabilizing effects of parasites on the underlying consumer-resource population cycles inherent to the herbivore-autotroph submodel (i.e., the classic Rosenzweig-MacArthur predator-prey model). Second, we assessed the extent to which parasite virulence (i.e., the degree to which parasites harm the host) influences the magnitude of trophic cascades in the caribou-Ostertagia-producer system by varying the negative effects of the parasite on host survival (α S ), fecundity (α F ), and resource intake (α I ). Last, we conducted a sensitivity analysis of all model parameters to assess the robustness of our inferences regarding parasite harm parameters and the general relationship between all model parameters and both the stability of the model resource-host-parasite system and trophic cascade strength.
Effects of parasitism on host and resource dynamics are host trait dependent. In the absence of parasites, the model collapsed to the Rosenzweig-MacArthur model for producers and herbivores, producing stable coexistence at low producer growth rates and low herbivory rates but oscillating herbivore population sizes as producer productivity increased (63). Therefore, as expected, parasite-free simulations using parameter estimates for the Ostertagia-caribou system displayed limit cycles ( Fig.  2A). Adding a parasite with only intensity-dependent negative effects on host survivorship stabilized the producer-herbivore interaction ( Fig. 2B; α S = 10 À3 , α F = 0, α I = 0), consistent with previous results from similar models (36). Adding a parasite with only intensity-dependent effects on host resource intake also stabilized the underlying producer-herbivore interaction ( Fig. 2C; α I = 10 À3 , α S = 0, α F = 0). Importantly, the stabilizing effects of parasite-induced declines in host resource intake occurred via a different mechanism than that of parasite effects on host survivorship (Fig. 2B). In this case, large parasite populations suppressed further parasite transmission by reducing host resource intake, which consequently reduced exposure of hosts to parasites. Finally, adding a parasite with only intensity-dependent effects on fecundity failed to stabilize producer-herbivore cycles, thus enabling the parasite population to grow exponentially ( Fig. 2D; Our model also demonstrates that negative effects of parasites on host traits can have qualitatively different dynamical consequences depending on the trait involved. Specifically, negative effects on host survivorship caused density-dependent increases in parasite mortality, because parasites within infected hosts are simultaneously killed with the hosts. Likewise, negative effects on host resource intake caused density-dependent decreases in parasite transmission. In both of these scenarios, stabilization occurs due to negative density dependence on the parasite. When mean infection intensity is high, average parasite mortality becomes higher and/or average exposure rate becomes lower due to the infection negatively influencing host survival or ingestion rate, respectively. In contrast, negative effects on host fecundity did not impose any density-dependent negative feedback on the parasite population and, therefore, did not stabilize the underlying producer-herbivore cycles. Interestingly, in simulations with parasite-induced negative effects on two host traits simultaneously, we found that effects on fecundity always destabilized producer-host-parasite dynamics, while negative effects of parasites on host survival or host resource intake consistently stabilized producer-host-parasite dynamics (Fig. 3). When parasite effects on fecundity were paired with one of the stabilizing effects (survival or resource intake), greater reductions in fecundity were required to generate cycles ( Fig. 3 A and C).
Our global sensitivity analysis indicated that these results generally hold over a much broader range of parameter space (Latin hypercube sampling sensitivity analysis on the presence and amplitude of cycles (64); ranges: 0 to 1 for all parameters representing proportions or probabilities; 10 À8 to 10 À3 for parasite harm parameters α I , α S , and α F ; 10-fold variation above and below estimates indicated in Table 1 for all other parameters). In addition, each mechanism of parasite-induced harm exhibited a different pattern (SI Appendix, Fig. S6). Parasite harm to resource intake was identified as consistently and strongly stabilizing to the producer-host-parasite system. In contrast, parasite harm to fecundity was consistently destabilizing. Finally, harm to host survival was stabilizing when parasites were aggregated [as is consistently seen for macroparasites (65)] but destabilizing when parasite aggregation was weak. Reductions in parasite aggregation among hosts also tended to destabilize the producer-host-parasite system (SI Appendix, Fig. S6).
Trophic cascades can be triggered by lethal and sublethal infections. Trophic cascades were triggered by lethal (survival) and both reproductive (fecundity) and nonreproductive (resource intake) sublethal effects of parasitic infection in our model. Specifically, in comparison to the parasite-free scenario (Fig.  4A), producer biomass was higher when parasites were present regardless of the host traits affected. Similar to classic predatorinduced trophic cascades, lethal parasitic infections reduced herbivore densities and indirectly released producers from topdown regulation (Fig. 4B). In addition, sublethal effects that reduced rates of host resource intake (Fig. 4C) or fecundity (Fig. 4D) increased producer biomass compared to the parasitefree scenario. Interestingly, trophic cascades triggered by sublethal effects of parasites were as strong as those triggered by lethal effects. One potential explanation of this phenomenon is that, all else equal, high virulence on host survival limits the abundance of parasites, thereby limiting their population-level impacts. An identical phenomenon can be seen in the R 0 criterion from classic susceptible-infected-type models, in which high host mortality limits the infectious period and thus total transmission potential (66).
In our global sensitivity analysis, we found that the magnitude of the trophic cascade was insensitive to the parasite harm parameters but highly sensitive to the foraging traits of the hosts (indicated by the relative change in plant biomass with vs. without the parasite; SI Appendix, Fig. S6). Specifically, larger trophic cascades occurred when hosts had higher feeding rates (f) and fecundity conversion efficiency (y) as well as lower background death rate (d H ) and lower half-saturation constant in their type II functional response (h A ). Taken together, these traits all contributed to a greater ability of hosts to suppress producer density in the absence of parasitism. Thus, larger trophic cascades can be triggered by parasites causing any type of harm if they release producers from regulation by rapidly feeding and efficient herbivore hosts [i.e., those with low R* values in the context of resource competition theory (67)].
Meta-Analysis Results. We used a meta-analysis to quantify the effects of helminths on the same herbivore traits (survival, fecundity, and resource intake) in natural systems as those shown to trigger trophic cascades in our model. To do this, we conducted a systematic literature search for studies that measured the effects of helminths on these three traits, in addition to body mass and body condition, in free-living and/or wild ruminants. After excluding studies which did not meet our inclusion criteria (SI Appendix, Fig. S1), our final dataset included 259 records from 59 studies spanning 18 host species (SI Appendix, Fig. S2) and five global regions (SI Appendix, Fig. S3).
Helminth effects on free-living ruminant hosts tend to be sublethal. We observed significant heterogeneity across all effect sizes (I 2 = 0.97, Q 258 = 14,134, P < 0.001), defined as the standard mean difference in each trait between infected and uninfected hosts. However, on average, helminth infection significantly reduced host performance as defined by the focal host traits (d = -0.47, 95% CI = -0.77 to -0.18, P = 0.001), especially for experimental studies (d = -0.65) compared to observational studies (d = -0.38; Q 1 = 3.89, P = 0.05, R 2 = 0.03). Examining each host trait individually, we found that helminth infection significantly reduced body mass (d = -0.61, 95% CI = -1.11 to -0.14) and body condition (d = -0.34, 95% CI = -0.66 to -0.03), and we found weaker but still generally negative effects of infection on feeding rate (d = -0.48, 95% CI = -1.03 to 0.03; Fig. 5 and SI Appendix, Table S2). Helminth infection was not significantly associated with either fecundity (d = -0.08, 95% CI = -0.21 to 0.05) or survival (d = -0.72, 95% CI = -2.05 to 0.40), although estimated mean effect sizes were both negative, and effect sizes for survival were highly variable (I 2 = 0.99). Additionally, across all host traits, the majority of records had effect sizes that were negative in direction (mass: 78%, body condition: 66%, feeding rate: 95%, fecundity: 82%, survival: 62%). The direction of these mean effect size estimates was generally consistent when we stratified this analysis by study type (experimental vs. observational), although sample sizes were small for some comparisons. When considering experimental data alone (SI Appendix, Fig.  S4 and Table S2), infection significantly reduced body mass and body condition; both feeding rate and fecundity were weakly negatively affected by infection as well. When considering just observational data (SI Appendix, Fig. S4 and Table S2), only fecundity and body mass were significantly and negatively associated with helminth infection. For those records that included caribou as the focal host (n = 71 records from 10 studies; SI Appendix, Fig. S5 and Table S2), helminth infection was negatively associated with body condition (d = -0.23, 95% CI = -0.43 to -0.03) and feeding rate (d = -0.87, 95% CI = -1.10 to -0.64), supporting the use of the caribou system as a case study for our modeling exercise.

Discussion
By combining a mechanistic model and meta-analysis, we found that both lethal and sublethal effects of helminth parasitism can modify interactions between herbivores and their resources but that sublethal effects occur more commonly in nature. First, our model showed that parasites that affect host survival or feeding rate stabilized producer-herbivore population cycles (Figs. 2 B and C and 3B), whereas parasites that affect host fecundity were more likely to have destabilizing effects (Figs. 2D and 3 A and C). Second, similar to predators, lethal effects of parasites (i.e., reductions in survival) triggered classic trophic cascades (2, 68) via a reduction in herbivore density and concomitant release of producers from top-down control (Fig. 4B). Third, sublethal effects of infection-including reduced host feeding rates (Fig. 4C) and reduced fecundity (Fig. 4D)-also caused trophic cascades, with effects on primary producer biomass that were as strong as those caused by reductions in herbivore host survival (Fig. 4B). Whereas parasites that reduced herbivore host resource intake triggered trophic cascades through effects on a host trait (i.e., feeding rate), parasites that affected host fecundity and survival triggered cascades through effects on host density (i.e., by reducing the number of herbivores on the landscape). Finally, while our model revealed that helminth parasites can trigger trophic cascades via diverse mechanisms, the results of our metaanalysis indicated that the actual effects of helminth infections on free-living ruminant herbivores are predominantly sublethal (Fig.  5). Given that free-living ruminants play keystone roles in many ecosystems (69)(70)(71) and that helminth parasites are ubiquitous, these sublethal effects (e.g., reduced feeding rates) could have broad-scale ecosystem consequences. Taken together, our findings support a role for sublethal parasitic infections as a global driver of ecosystem dynamics through the modification of host functional traits that ultimately impact producer biomass. Interestingly, when varying the intensity of multiple parasite effects on hosts simultaneously, we found that certain modes of parasite harm can stabilize host-resource oscillations, whereas others may destabilize host-resource interactions ( Fig. 3 and SI  Appendix, Fig. S6). The stabilizing effects of harm to host feeding rate that we observed for the Ostertagia-caribou system stand in contrast to the destabilizing effects of harm to fecundity both in this system and in models of helminth infections in red grouse populations (72). Our model indicates that these alternative stabilizing and destabilizing outcomes arise via different density-dependent mechanisms. Specifically, while The model was parameterized with values corresponding to the producer-caribou-Ostertagia system (Table 1) with herbivores, adult parasites, and freeliving stages of the parasite expressed as number per square kilometer and producer biomass expressed in kilograms per square kilometer. The density values are visualized as relative to the maximum of each group in each panel.
negative effects of parasites on host fecundity do not directly impose any density-dependent negative feedback on the parasite population, infection-induced reduction of host feeding rate does cause density-dependent decreases in parasite transmission because hosts are exposed to parasites through feeding. We also found parasite harm to host survival can be stabilizing, but only when parasites are highly aggregated. High levels of parasite aggregation have similarly been found to promote stability in models of red grouse-nematode dynamics (72). Aggregation is typical for macroparasites (65), including for Ostertagia in caribou (73,74). Consequently, a mechanistic understanding of how parasite aggregation interacts with different modes of harm to host traits is critical for predicting the stability of different resource-host-parasite systems. Notably, irrespective of whether effects of parasites on host-resource interactions were stabilizing or destabilizing, trophic cascades were possible. Fig. 3. Two-parameter diagrams illustrating the qualitative behavior of the mechanistic model (Fig. 1), which oscillates in the absence of disease. Panels display variation in the magnitude of parasite virulence for pairs of lethal (survival) and sublethal (fecundity, resource intake) components included in the model (log 10 -transformed, setting the third component to zero). Parasites that harm any combination of host survival and resource intake consistently stabilize oscillations (gray; B). In contrast, parasites that harm host fecundity generate producer-host-parasite oscillations (orange), even in combination with a component of virulence that is otherwise stabilizing (A and C). The model is parameterized with values corresponding to the producer-caribou-Ostertagia system (Table 1)  Indeed, our mechanistic model indicates that trophic cascades caused by helminths may occur much more commonly than previously recognized. Specifically, both lethal and sublethal effects of infection readily generated trophic cascades across a gradient from low to high virulence. For parasites that reduced host resource intake, the strength of the trophic cascade was consistent across the virulence gradient. However, increased negative effects of parasites on host survival resulted in a very slightly weakened trophic cascade on producer biomass (not detectable in Fig. 4B) due to an increase in host density. This pattern of increased host density occurred in part because higher host mortality shortened the infectious period, resulting in reduced parasite population densities (75). Overall, however, we found that there was little sensitivity to the magnitude of virulence for any of the three host traits investigated. By contrast, in another recent modeling study, equilibrium biomass of plant hosts declined with increasing virulence of a pathogen (30). Fungal parasites of zooplankton hosts can also trigger trophic cascades or lead to increased host density, depending on resource productivity and the relative harm of infection to host feeding rate versus survival (29). These divergent outcomes highlight the need for future work delineating how parasite effects on ecosystems vary according to the type of parasite (e.g., macro-vs. microparasites), host trophic level, and type of harm to hosts. Interestingly, our results contrast with previous work suggesting that for macroparasites aggregation reduces their potential to trigger trophic cascades (2). Instead, we find that these common parasites have broad potential to induce ecosystem-level effects.
Despite the potential for helminths to alter community dynamics through multiple mechanisms, nonreproductive sublethal effects appear to be the most common way in which helminths trigger trophic cascades through ruminant hosts in nature. Our meta-analysis indicated that helminth infection consistently reduced host feeding rate, body condition, and body mass but on average did not affect fecundity or survival. Persistent negative effects of helminths on host body condition and mass may reflect reduced resource intake by infected individuals, but downstream effects on demographic traits such as host fecundity or survival may be weaker or more difficult to detect. Overall, however, these results support the expectation that helminth infections are harmful to host health but that effects detected in free-living ruminants are predominantly sublethal (76). To date, studies on parasite-mediated trophic cascades have focused on the lethal effects of viruses and other microparasites and the sublethal, reproductive effects of macroparasitic castrators (2,3,7). Our findings show that noncastrating macroparasites can also trigger trophic cascades through nonreproductive consumptive effects on hosts (2). Given that nonreproductive sublethal consequences of helminth infection are common in many host-parasite systems, our findings highlight the need for further research on the community-and ecosystem-level consequences of diverse parasite effects on host function. In particular, parasite-induced changes in host feeding behavior, which have been documented across a range of host and parasite taxa (60), could have widespread ecosystem impacts. It is important to note, however, that although helminth infection was not associated with fecundity or survival on average in our meta-analysis, variation in effect sizes was extremely high for survival; in fact, several studies showed negative effects of parasitism on survival. The majority of all effect sizes per trait also showed negative signs regardless of statistical significance. Our mechanistic model suggests that even small negative effects of infection on feeding, fecundity, or survival can have important consequences for host population dynamics and producer biomass. Thus, in contexts where helminth infection reduces host survival or fecundity by even a small amount, as observed in a large fraction of studies, there may still be cascading effects on ecosystems via shifts in producer biomass production.
The findings from this study suggest that there are global impacts of sublethal parasitic infections on ecosystems. Large herbivores often play an outsized role in ecosystem function (69,70,77), especially in ecosystems with few trophic levels (54,78,79); thus, even subtle parasite impacts on herbivore traits may have cumulative cascading impacts on ecosystems. Given that helminth parasites are nearly ubiquitous within free-living populations of ruminants (44,46,80), and that our meta-analysis of helminth parasitism in ruminants was global in its coverage, we suggest that global herbivory rates by ruminants are likely lower as a consequence of pervasive helminth infection. If so, the world may be greener due to the top-down effects of parasites on herbivore hosts. For example, for studies in our meta-analysis in which caribou were the host species, helminth infections caused significant reductions in host feeding rates, body mass, and body condition (SI Appendix, Fig. S5). Furthermore, our model showed that producer biomass increased when caribou were infected with Ostertagia, a common gastrointestinal parasite in this species (48).  Thus, in systems like the Arctic, parasites may play a significant and underappreciated role in shaping producer dynamics. Experimental work linking hosts, parasites, and producers will be required to understand if the cascading effects of parasites on ecosystems are similar in magnitude to those induced by some predator-prey interactions (e.g., refs. 1 and 81). Future work would benefit from the development of stoichiometrically explicit models [e.g., as recently developed for plant pathogens by (30) and predator-prey systems by (82)] and collection of relevant data to clarify how sublethal effects of parasites on herbivorous hosts alter fluxes of carbon and elemental nutrients through ecosystems and how such impacts vary across ecosystems and biomes.

Conclusion
The potential for parasites to impact ecosystems, and vice versa, has been a topic of interest in recent efforts to link the fields of disease ecology and ecosystem ecology (2,3,7,30,38,83,84). One of the major limitations in predicting parasite effects on ecosystems is that while studies of how parasites influence host individuals or populations have been common, syntheses of parasite effects on specific host traits that can be linked to ecosystem-level processes are rare. By using a mechanistic model to connect individual-level effects of infection to population and community dynamics and a meta-analysis to evaluate empirical support for these trait-based patterns in nature, our approach furthers our understanding of the epidemiological and ecosystem consequences of common parasitic infections. Characterizing the broader ecological roles of parasites is increasingly important, because their effects on hosts in both natural and agricultural systems are in flux. For example, parasite abundance among some of the world's most common livestock animals is increasing (85,86), while similar parasites in wildlife are facing extinction within the next half century (9). Our findings reveal that such changes in parasite abundance and distributions may have widespread cascading impacts on ecosystems.

Materials and Methods
Mechanistic Model. We tested the potential of parasites to trigger trophic cascades by using a mechanistic model of macroparasite (helminth) infections in an herbivorous host that builds on previous work (36). We used ordinary differential equations to track the densities of autotrophic producers (A), host/herbivores (H), adult parasites within hosts (P), and free-living stages of the parasite in the environment (E) through time (Fig. 1). We illustrated several key features of the model and obtained general results by parameterizing the model with data from a plant/lichen-caribou-parasite (Ostertagia) system in the arctic tundra (Table 1).
Producers grow logistically, with maximum growth rate, r, and carrying capacity, K, reflecting competition for light or nutrients. Herbivores consume producers with a type-II functional response with maximum foraging rate, f, and half saturation constant, h A . Parasites can exert intensity-dependent negative effects on host resource intake at cost α I (Eq. 1). Herbivorous hosts reproduce at a total rate that depends on ingested resources multiplied by their fecundity yield, y, such that reduced resource intake due to parasitism has a downstream negative effect on net fecundity of hosts. Parasites can also harm fecundity yield with intensity-dependent effects at cost α F . Hosts die with a background death rate, d H , and also due to intensity-dependent negative effects of the parasite on survival, α S (Eq. 2). Adult parasites in hosts, P, increase due to infection, which occurs with per-parasite infection probability, σ, following ingestion of free-living parasite propagules (inclusive of eggs and larvae), E, and is linked to the host's type-II functional response to producers. Adult parasites are lost due to their own background mortality, d P , and when hosts die naturally or due to virulent effects on survival (Eq. 3). To avoid the possibility of negative rates of resource intake or biomass yields, which are biologically impossible, we represent these effects as relative decreases in these trait values per parasite; thus, the α F and α I terms represent the relative decrease in these two processes as infection intensity increases. We followed ref. 87 to integrate this assumption of relative decrease with the negative binomial distribution of parasites among hosts to generate terms representing the total ingestion and fecundity across the host population (terms raised to the exponent k in Eqs. 1-4). The final term in Eq. 3 represents the effect of parasite aggregation on the loss of parasites when heavily parasitized hosts die. Parasite propagules in the environment are produced at rate ρ by adult parasites, die at background rate d E , and are removed via ingestion by hosts during the transmission process. We explored the effects of parasites on the dynamics of this producer-herbivore system by solving the model numerically using the deSolve package in R (88,89). We explored variation in each parasite-induced effect on hosts (survival [α S ], resource intake [α I ], and fecundity [α F ]), first individually (while setting the others to zero) and then in all pairwise combinations. We ran each simulation for 40,000 daily time steps (∼110 y); after allowing for initial transient dynamics, we evaluated long-term dynamics (i.e., stable equilibrium values or the mean and SE of each state variable if the system cycled) using the last 10,000 time steps. We conducted a global sensitivity analysis of the model to further evaluate the generality of our inferences regarding parasite effects and stability. First, we focused on two endpoints, characterizing either the presence of cycling (0 = no, 1 = yes) or the magnitude of cycling, using the coefficient of variation of host density. For the third endpoint, we focused on the magnitude of the parasiteinduced trophic cascade, with the endpoint represented as proportional change in plant abundance with and without the parasite. We used a Latin hypercube sampling design with the randomLHS() function from the lhs R package (64). We drew 60,000 randomly generated parameter sets across broad ranges of parameters and then simulated the model as above. Parameter ranges for this sensitivity analysis were 0 to 1 for all parameters representing proportions or probabilities and 10 À8 to 10 À3 for parasite harm parameters α I , α S , and α F . All other parameters varied 10-fold above and below the estimates indicated in Table 1. We sampled the parasite harm parameters uniformly over logtransformed space to more evenly sample the 10 5 -fold gradients. We then excluded any parameter sets in which the parasites failed to invade the population (∼12% of random parameter sets), because in these cases parasite traits would be irrelevant. We then computed partial rank correlation coefficients (PRCC) for all model parameters with respect to both measures of cycling. PRCC values span from À1 to 1 and can characterize monotonic relationships between model input parameters and model outputs without assuming linearity (90). In this case, larger values of either model output indicate cycling (destabilization), and therefore a negative PRCC value for a parameter represents evidence that increases in that parameter contribute to stabilizing the resource-host-parasite dynamics whereas a positive PRCC value indicates that increasing that parameter destabilizes the system. Datasets S5 and S6 contain the R code for the mechanistic model and sensitivity analysis.

Meta-Analysis.
Literature search. We performed systematic searches in BIOSIS, Science Citation Index, Zoological Record, and PubMed for papers published between 1 January 1980 and 24 February 2019 (the date of searches) using search strings restricted to free-living ruminant-helminth interactions (SI Appendix). We retrieved a total of 4,371 records, 2,217 of which were unique records after removing duplicates based on matching title, abstract, and/or DOI. We screened the titles and abstracts of all 2,217 unique records using metagear (91) in R and excluded studies that did not pertain to free-living ruminant-helminth interactions (i.e., all domestic or livestock hosts were excluded) or that did not include data on any of the focal host variables (i.e., survival, fecundity, feeding rate, body mass, or condition). We retained 390 studies after this first screening phase and assessed the full text of each of these articles for eligibility. In several cases where we determined that a study was potentially eligible but lacked necessary information for inclusion, we contacted the authors directly to request clarification. Ultimately, 59 articles contained sufficient data and methodological details for inclusion in the meta-analysis (SI Appendix and SI Appendix, Fig. S1). Data extraction. From each study, we recorded the ruminant species, helminth species (or broader taxon), response variable (i.e., body mass, body condition, fecundity, feeding rate, survival), sample size, test statistics for the relationship between parasitism and the response, and the direction of effect (i.e., positive or negative) between parasitism and the response. We also recorded study attributes including location, type (i.e., experimental or observational), and duration.
We reported effect size as the correlation-based r between parasitism and host traits (92). Many studies did not directly report correlation coefficients, so we first converted other test statistics (e.g., χ 2 , F, t, z, odds ratios) into r (93,94). When test statistics and effect directionality were not reported, we derived them using summary data (e.g., calculating Cohen's d based on means and SDs or odds ratios from tabulated counts) or converted the P value to a standard normal deviate Z-score and used the sample size to obtain r. When SDs were not reported, we estimated these using ranges or 95% confidence intervals (95,96). When dispersion measures were unavailable, we imputed missing SDs using the weighted mean per study or across all data (when a study was missing one variance estimate). When exact P values were not reported (e.g., P < 0.05), we followed guidelines to convert p to Z (92). We assigned negative values to r when host traits were lower in infected animals. Using the metafor package, we converted all directional r into Fisher's Z (Z r ) (97). Records were excluded if these data were not reported or could not be derived or if sample sizes were less than four. The data for all records included in the meta-analysis are available in SI Appendix. Meta-analysis. To assess the overall relationship between helminth infection and ruminant traits, we used hierarchical meta-analysis models with the metafor package (97,98). We nested observations within a study-level random effect to account for unit-level variance and pseudoreplication, as most studies had multiple effect sizes (34/59). We compared a series of additional random effects in order to account for repeated measures in a subset of studies (6/59), phylogenetic dependence of host species (studies did not consistently report helminth species to allow a random effect for parasite phylogeny), and spatial autocorrelation, with restricted maximum likelihood (99)(100)(101)(102). A structure with only observation, study, and an autoregressive structure was the most parsimonious, although all possible random effects generated similar estimates of model coefficients (SI Appendix, Table S1). In particular, we did not find additional support for phylogenetic random effects, which suggests ruminant phylogeny does not have strong impacts on effect size and that relationships between helminths and host traits may be relatively homogenous across ruminants. All models used the rma.mv () function, weighting by sampling variance, and were fit with restricted maximum likelihood (REML). In all analyses, we back-transformed Z r to r and then to d for interpretability (i.e., the standardized mean trait difference between infected and uninfected hosts).
To assess heterogeneity among all effect sizes, we fit a random-effects model (REM; intercept only). We used the estimates of the variance components from this REM to quantify I 2 , the contribution of true heterogeneity to the total variance in effect size (103). We used Cochran's Q to test if such heterogeneity was greater than that expected by sampling error alone (93). Next, to assess whether effect sizes varied between experimental and observational studies (104), we fit a mixed-effects model with the same random effects and study type as a moderator. We assessed moderator significance using the Q test (97) and derived a pseudoR 2 as the proportional reduction in the summed variance components compared with those of the equivalent REM (105).
To estimate the mean relationship between helminth infection and each host trait, we repeated the above analysis by fitting a REM with the same random effects structure to each host trait subset of our data. Because the sample size for several of these traits was relatively low, we a priori stratified these analyses by observational and experimental study to assess whether stronger effects of parasitism on each trait were detected using experimental approaches (104).